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High- frequency power amplifiers operate most efficiently at satu- 
ration, i.e., in the nonlinear range of their input /output characteris- 
tics. This phenomenon has traditionally dictated the use of constant 
envelope modulation methods for data transmission, resulting in 
circular signal constellations. This approach has inherently limited 
the admissible data rates in digital radio. In this paper we present a 
method for solving this problem without sacrificing amplifier power 
efficiency. We describe and analyze an adaptive linearizer that can 
automatically compensate for amplifier nonlinearity and thus make 
it possible to transmit multilevel quadrature amplitude modulated 
signals without incurring intolerable constellation distortions. The 
linearizer utilizes a real-time, data-directed, recursive algorithm for 
predistorting the signal constellation. Our analysis and computer 
simulations indicate that the algorithm is robust and converges 
rapidly from a blind start. Furthermore, the signal constellation and 
the average transmitted power can both be changed through software. 



I. INTRODUCTION 

Progress in high-speed data transmission over radio channels has 
lagged behind that of the voiceband channel. Inherent difficulties 
associated with implementing automatic equalizers are partially re- 
sponsible, but the application of multilevel quadrature amplitude 
modulation (QAM) also has been inhibited by the amplitude 
(AM/AM) and phase (AM/PM) nonlinearities present in radio-fre- 
quency (RF ) power amplifiers. 

Recent work 1 has evolved design principles showing a possibility of 
substantial improvements in QAM performance over linear fading 
radio channels. A crucial obstacle to achieving these gains is the 
nonlinear distortion introduced by power amplifiers. Attempts to re- 
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alize high-speed data transmission over these channels force consid- 
eration of methods to cope with this nonlinear distortion. 

One approach is to back off from saturation sufficiently so that the 
signal level is restricted to the linear range of amplification. The 
required amount of power back-off can be several decibels, resulting in 
an inefficient operation of the power amplifier. Moreover, the achieve- 
ment of a given desired level of average transmitted power would 
require the use of a large, expensive, and high-power-consuming am- 
plifier. 

It has been realized 2-7 that some improvements in this regard can be 
obtained by using fixed signal predistortion circuits prior to amplifi- 
cation. Such circuits, however, cannot compensate for drifts in power 
amplifier nonlinearities caused by temperature changes, dc power 
variations, and component aging. These fluctuations can considerably 
degrade performance of systems employing constellations with large 
numbers of points, say 64 or higher, and the use of an adaptive 
technique is necessary in these applications. 

Conceptually, the nonlinear distortion introduced by the power 
amplifier can be minimized at the receiver by an adaptive nonlinear 
equalizer. Such schemes have been proposed and studied in voiceband 
data transmission 8 and for filtered PSK signals operating over satellite 
channels. 9 This approach does not seem to be reasonable or necessary 
in our application, since the source of the nonlinearity is at the 
transmitter. Thus, it would appear logical to equalize the nonlinearity 
at the transmitter, where it occurs and where the transmitted bits are 
available. 

Thus, this paper focuses on the problem of adaptive predistortion 
linearization. We describe a transmitter-based recursive algorithm for 
predistorting the signal constellation, thereby rendering a virtually 
linear transmitter. The algorithm operates in real time and is data 
directed. The predistortion is accomplished within a digital memory, 
which is used to generate the desired baseband signal. This maximizes 
the use of digital technology, and increases the reliability and flexibility 
of the system. For example, it is possible to change the signal constel- 
lation and the average transmitted power through software. Our treat- 
ment applies only to single-valued, memoryless nonlinearities. 

The idea of adaptive predistortion of signal constellations has been 
previously suggested. 10 In this reference, the predistortion of each 
point of the constellation is accomplished by switching the RF or 
intermediate-frequency (IF) input signal to a separate path containing 
an adjustable diode attenuator and an adjustable diode phase shifter. 
The analog hardware required for such an implementation would be 
quite involved, especially for a large number of points. We have also 
found in the patent literature 11 a description of a digital adaptive 
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predistorter that appears to be similar to the one we present here. 
However, no detailed description of the algorithm or its behavior is 
provided. 

Our basic ideas are described in Section II. Mathematical analysis 
is provided in Section III, and simulation results of the adaptive 
linearizer with finite-precision arithmetic are given in Section IV. 

II. GENERAL DESCRIPTION AND REQUIREMENTS OF OPERATION 

A block diagram of a QAM transmitter with the proposed predistor- 
tion linearizer is shown in Fig. 1. A random access memory (RAM) 
contains the predistorted values of the in-phase and quadrature volt- 
ages of each point on the QAM constellation. A memory-lookup 
encoder obtains each input data symbol and generates the RAM 
addresses of the desired signal point. The corresponding stored, pre- 
distorted voltage values are converted to analog voltages using a pair 
of digital-to-analog (D/A) converters. These voltages drive a quadra- 
ture modulator, which generates the desired predistorted RF signal for 
the duration of the input symbol. That signal is then amplified, filtered, 
and transmitted. This part of the linearizer is similar in operation to a 
recently proposed memory-based encoder described in Ref. 12. 

To update the RAM information, the amplifier output is sampled 
by a directional coupler and demodulated using the same local oscil- 
lator used for the modulator, which eliminates the need for carrier 
recovery. The output in-phase and quadrature voltages of the de- 
modulator are converted to digital form using a pair of analog-to- 
digital (A/D) converters. A linearizing processor, which is the heart of 
the linearizer, receives this information, compares it with the input 
data, and computes the resulting error. A recursive algorithm, which 
is discussed in the next section, uses the error to update the voltage 
values in the RAM corresponding to the particular data point under 
consideration. Note that each point on the signal constellation is 
treated separately. Thus, the linearizer can support any desired con- 
stellation. 

The memory-lookup encoder and the D/A converters have to op- 
erate, of course, at the full signaling rate. However, the linearizing 
processor and the A/D converters can operate at a much reduced rate, 
since the updating process is only needed to compensate for drifts that 
occur on a much slower time scale than the data rate. 

We now emphasize a crucial point. The proper operation of the 
linearizer as described above assumes that the amplifier is memoryless, 
and requires that the signal not be filtered before the power amplifier. 
Thus, all pulse shaping and filtering must be performed by the com- 
bination of the post- amplifier, RF bandpass filter, and the receiver 
filters. The former filter should be designed just to meet FCC (or 
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other) emission rules for square input pulses, as described in Ref. 13. 
Such an implementation may require an automatic equalizer at the 
receiver to eliminate residual intersymbol interference. It was shown 1 
that, even with ideal filtering, adaptive equalization would still be 
necessary to compensate for multipath fading using QAM signals with 
large numbers of levels. Thus, the elimination of pre-amplifier pulse 
shaping appears to be a mild constraint. 

Because of the more stringent requirement on the post-amplifier, 
RF bandpass filter, its loss would of course be increased over that of 
the conventional case where preamplifier filtering does the spectrum 
shaping. However, computations 14 based on practical filters 15 operating 
in the 6-GHz band show that the RF filter loss would only increase 
from about 0.5 dB to 1.5 dB. However, we will see in Section IV that 
the linearizer in our system would allow the operation of the power 
amplifier to approach saturation. This would result in several decibels 
of power increase, which would more than compensate for the one- 
decibel increase in the filter loss. 

III. THE RECURSIVE ALGORITHM 

As we already mentioned, our approach is applicable to any signal 
constellation. However, for clarity of exposition, we restrict our anal- 
ysis to a rectangular constellation. Referring to Figs. 1 and 2, we denote 
a point on the rectangular QAM constellation by the complex number 

a + ib = pe*, 

where, for an L 2 -level system, a and b assume values on the lattice 
±1±3±5± ••• ± (L - 1). We denote the predistorted point by the 







" 


• 


• 


® • 


• 


• 


• • 


• 


• 


• • 


• 


• 


• • 



♦ ^ 



® 



® 



QUADRATURE AMPLITUDE 


DIGITAL/ANALOG 


ANALOG/DIGITAL 


MODULATION DATA 


DATA 


DATA 


a + ib = pe'* 


x + iy = re'" 


u + ivs Re'* 



(a) 



(b) 



(c) 



Fig. 2 — The signal constellations at various points in the linearizer of Fig. 1. 
(a) Quadrature amplitude modulation data, (b) Digital/analog data, (c) Analog/digital 
data. 
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complex number 



x + iy = re iB . 



A sequence of these points amplitude and phase modulates a carrier at 
frequency f c and the resulting signal, 



S(t) = £ r< m) p(t - mT)e i[2,rf < t+6{m)] , 



(1) 



is applied to the amplifier. In (l),p(t) is a rectangular pulse, 1/Tis the 
signaling rate, and r (m) , (m) represent the amplitude and phase of the 
with data symbol. 

The amplifier is customarily represented by a pair of memoryless 
nonlinear functions. 16,17 The amplitude function, A{r), causes ampli- 
tude distortion (AM/ AM) and the phase function, O(r), causes ampli- 
tude to phase conversion (AM/PM). Thus the signal, (1), after ampli- 
fication becomes 



S (t) = £ A(r {m) )p(t - m T)e il2 < t+eim)+ ° lr(m)] \ 



(2) 



Figure 3 shows sketches of typical curves A (r) and <I>(r) for a traveling- 
wave tube (TWT) power amplifier. 

We remark that if the functions A ( • ) and <D( • ) were known exactly 
we could choose transformations g(>) and h(-) from points (p, <>) to 
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Fig. 3 — Typical example of the input/output amplitude and phase characteristics of 
a TWT amplifier. 
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(r, d) so that 



where 



and 



where 



A(r) = A[g(p)] = Gp, (3) 

r = g(p) 
$(r) + h(<f>, p) = <t> + f, (4) 

= n(4>, p). 



The constant gain, G, in (3) is some desired gain, taking into account 
the linear gain of the power amplifier and the coupling factor of the 
sampling directional coupler. The fixed phase shift, f , in (4) is arbitrary 
and can be set to zero without loss of generality. 

Finally, let the measured data point from the A/D converter be 

u + iv = Re 1 *. 

The object of the predistortion is to find a solution 

r = f 

and 

= 6, 

such that 

Re 1 * = Gpe*. (5) 

We now describe an iteration procedure that converges to (5). 

Since data are usually scrambled, a specific data point (p, <>) will 
occur at random. Let (r„ , 6 n ) be the predistorted RAM data point and 
(R n , »M be the measured data point at the nth time the desired data 
point (p, <f>) occurs. The measured output radius R n is then 

Rn = A{r n ) + v n (6a) 

while the measured angle is 

^ n = 0(r„) + B n + tin, (6b) 

where v n and ju n are zero-mean measurement errors. If the measure- 
ments were perfect and noiseless, we would simply solve the following 
two nonlinear equations 

f(r)=A(r)-Gp = (7) 
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and 

$(r) + 0-4> = O (8) 

for r and 9. 

A great number of iterative procedures are known for solving the 
nonlinear equations given in (7) and (8). In the presence of measure- 
ment noise (with convex functions A and O), stochastic approximation 
algorithms 18 provide efficient methods. So, if one chooses step sizes 
a n and /?„, which behave as o(l/n), the following recursions are known 
to converge to the true solutions f and 6 in a mean-square sense: 

r„ + i = r„ - a n (R n - Gp) (9a) 

e n+ i = 9 n - fintyn - <f>). (9b) 

To provide a rationale for the above recursions we analyze, in some 
detail, the behavior of (9a). The behavior of (9b) is similar. Since f(X), 
eq. (7), is continuous on (Xo, Xmax) and if, by hypothesis, the derivative 
f'(X) of f(X) exists on this interval, then there is a £ such that 

Xo<£<Xma* 

m = f(X^)-f(Xo) m 

Amax -A0 

Applying this mean- value theorem to f(r„), eq. (7), we get 

f(r n ) = A(r n ) -Gp= f(f) + /'(^)(r„ - r), (11) 

where £„ lies between r„ and f. Substituting this into (9a) we get 

r„ +1 -r=r n -r- a„[A'(|„)(r„ - r) + v n \ (12) 

In (12) we made use of the fact that f(f) = and subtracted r from 
both sides of the equation. 
Letting a„A'(^ n ) = y„, r„ — r = S n and then iterating (12) we obtain 

On+l = (1 - yn)O n — "rrj- 



A'(tn) 



= 5. n (i - y.) - s t^t 1 - n (i - yi). d3) 

We now compute the mathematical expectation of 8 n and the variance: 
5„ +l = E{8 n+1 ] = EiSi) n (i - yd 

<£{5i}e~-> Y ' (14) 



and 



aiU, = E{8 n+i - 5 n+I } 2 = ol l ? - 2 f[ (1 - y,) 2 , (15) 

*_i [A \$k)\ (1 — yh) i=k 
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where o\ -* E{vl). We see from (14) and (15) that rapid convergence 
of the algorithm is critically dependent on the step-size sequence, y„, 
which in turn depends on the derivative of the nonlinearity in the 
neighborhood of the solution. While the derivative is unknown, in most 
applications it can be bounded away from zero and this makes it 
possible to estimate the best step-size sequence. We see from (14) that, 
even in the absence of noise, the algorithm is guaranteed to converge 
to the true solution only if the step-size series diverges, i.e., 

I y, , - £ atA'di) -* co, n -+ oo. (16) 

i-i /=i 

If a, and A'(£i) are restricted to be positive, (16) is equivalent to 
requiring £"-i a, to diverge. So it follows that the structure of the 
sequence can be of the form 

*n = -,, (17) 

n 
where a is a positive constant and < tj < 1. It can be shown that for 
this choice of a„, and for c = aA'(r), 



ii-i) 



Ei^) =e ~ + ^\ ( o<tj<1) (18) 

E{8r] 



and 



where 



.2 



[A'(r)] 2 ^ = JL (0<tj<1), (19) 



{n l ~ n " 1 1 
1 — Z~ ^ T* ; 
1 - TJ k -\ K 'J 

e.g., q(0) = and q(l/2) » 1.46. Table I shows the behavior of the 
statistics for the special cases 77 = and tj = 1, where (18) and (19) are 
not applicable. 



Table I — Behavior of statistics for tj = and tj = 1 
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It is noted from eqs. (17) through (19) and Table I that the choice of 
a fixed step size (tj = 0) gives the fastest convergence of the mean, 8 n +i , 
but results in a finite variance, a £+i , as n — > oo. On the other hand, a 
variable, progressively smaller step size (tj ^ 0) gives a slower conver- 
gence of the mean, but results in a variance approaching zero as n — > 
oo. Thus, a fixed step size should be used if the measurement error 
variance, al, is very small, and a variable step size should be used if 
the variance is large. Actually, when finite-precision arithmetic is 
employed, the choice of a progressively smaller step size can lead to 
large errors. 19 Thus, in our particular application where the measure- 
ment error can be made small, and where the use of finite-precision 
arithmetic is necessary, a fixed step size is more suitable. 

We note that our algorithm, (9), is based on the polar representation 
of the data symbols, while the hardware implementation of Fig. 1 is 
based on the rectangular representation. Thus, the linearizing proces- 
sor is required to convert back and forth between the two representa- 
tions. One could avoid this conversion by replacing (9) by the rectan- 
gular-based algorithm 

x n +i + iy n +i = x n + iy n - {a n + ip n )[u n + iv n - G(a + ib)]. 

Unfortunately, however, the convergence of this algorithm is not 
guaranteed, even for some well-behaved amplifier characteristics. 

IV. SIMULATIONS 

Here we present the results of computer simulations of the algo- 
rithm, (9), with fixed step sizes and finite-precision arithmetic. The 
precision is limited by the finite number of bits of the D/A and A/D 
converters (Fig. 1). As an example, we consider a TWT amplifier with 
normalized amplitude and phase nonlinearities of the form 17 

1 + r z 

•to = 60° -4-2, (20b) 

1 + r z 

which are sketched in Fig. 3. Note that at saturation, r = 1, A(r) = 1, 
and <D(r) = 30°. Figure 4 shows the severe distortion of the output 
constellation obtained with such an amplifier for a 64-QAM input 
signal driven with its corner point at saturation. 

Figures 5, 6, and 7 show the simulated results of the output signal 
constellations after the application of (9) for three different cases 
(explained below). In all cases the amplifier drive is maintained at 
saturation. The four quadrants in each figure represent different 
combinations of the numbers of bits of the D/A and A/D converters, 
as indicated in Table II. 

1028 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1 983 



Fig. 4 — The distortion of the signal constellation obtained with the amplifier of Fig. 
3 for a 64-QAM input signal with its corner points at saturation. 

Figure 5 corresponds to the case of no measurement error, i.e., v„ = 
/x„ = in (6). The step sizes used in (9) are a n = 0.5 and /?„ = 1.0. Since 
the normalized small-signal gain, A'{r), of the amplifier is about 2, the 
choice of a„ = 0.5 results in a value for c of about 1. This results in the 
fastest convergence of (9a), as can be seen from the first row of Table 
I. Similarly, the choice of ft n = 1.0 results in the fastest convergence of 
(9b) since that equation is linear in 6. The initial guesses of r and 6 for 
each constellation point in Fig. 5 were chosen at random, i.e., from a 
blind start. The results indicated in the figure are those after 25 
iterations for each constellation point. The point scatter shown is 
entirely due to the finite resolution of the D/A and A/D converters 
since no measurement errors were assumed. 

It is clear from Fig. 5 that the use of 8 bits for both the D/A and 
A/D converters (fourth quadrant) gives quite an acceptable perform- 
ance for 64 QAM. The use of 9 bits for each converter (first quadrant) 
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Fig. 5 — Simulated results of the output signal constellation after 25 iterations with no 
measurement errors, and for step sizes a„ = 0.5 and {!„ = 1.0. The number of bits of the 
D/A and A/D converters are different in each quadrant as indicated in Table 2. 

results in an almost perfect output constellation. In general, for L 2 - 
QAM, where L = 2* 1 , D/A and A/D converters with (M + 5) bits are 
needed for acceptable performance, and (M + 6) or more bits are 
needed for almost perfect performance. 

In Fig. 6, the step sizes employed are a„ = 0.5 and (S n = 1.0, as in Fig. 
5. However, a measurement error was introduced that is equivalent to 
a 30-dB signal-to-noise ratio at the corner points of the constellation. 
With the normalization used for A(r) in (20a), this noise corresponds 
to a 2 . = A 2 (r)al = 0.0005, where at and a 2 are the variances of the 
errors v n and fi n defined in (6). It is clear from Fig. 6 that such a level 
of noise, in combination with the large step sizes employed, gives 
unacceptable results. 

In Fig. 7, the same noise level as that of Fig. 6 is employed. However, 
the step sizes were reduced by a factor of 10, i.e., a n = 0.05 and fS n = 
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Fig. 6 — Same as Fig. 5, but with a measurement error equivalent to 30-dB of signal- 
to-noise ratio at the corner points of the constellation. 



0.1, resulting in a greatly improved performance. About 100 steps of 
iteration were needed in this case to reach convergence. Note that in 
spite of the measurement noise, the performance in the third quadrant 
of Fig. 7 is better than the corresponding performance in Fig. 5, where 
no noise is present. This is due to the reduced step sizes, and to the 
fact that the noise tends to smooth out quantization errors. 

V. SUMMARY AND CONCLUSIONS 

We have proposed and analyzed a transmitter-based, adaptive 
linearizer, which automatically compensates for power amplifier non- 
linearity in digital radio systems employing multilevel quadrature 
amplitude modulation. The linearizer utilizes a real-time, data-directed 
recursive algorithm for predistorting the signal constellation. The 
algorithm is robust and results in rapid convergence, even from a blind 
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Fig. 7 — Same as Fig. 6, but with the step sizes reduced by a factor of 10, i.e., a n = 0.05 
and fi n ■ 0.1. About 100 iteration steps were needed for convergence. 



Table II — Number of bits of the D/A 

and A/D converters for the four 

quadrants of Figs. 5, 6, and 7 



Number of Bits Number of Bits 
Quadrant of the D/A of the A/D 

Number Converter Converter 



start. The digital nature of the linearizer allows both the signal con- 
stellation and the average transmitted power to be changed through 
software. 

The proposed linearizer is only suitable for amplifiers with essen- 
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tially memoryless and single-valued nonlinearities. Consequently, 
pulse shaping is not permitted before amplification. In our approach 
the task of pulse shaping is relegated to the combination of the 
transmit RF filter, a receiver-based filter, and an automatic equalizer. 
The latter is presumed to be required anyway to deal with multipath 
fading, especially for signals having a large number of levels. 

Our principal conclusion is that high-power RF amplifiers can op- 
erate at saturation, where they are most efficient, provided that an 
adaptive predistortion linearizer is used. We have demonstrated that 
such a system is feasible and have provided a theory and a methodol- 
ogy for assessing its performance. 
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